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Abstract 

The existence of a shallow or virtual tetraquark state, ccud, is discussed. 
Using the putative masses for the doubly charmed baryons {ecu/ ccd) from 
SELEX, the mass of the ccud state is estimated to be about 3.9 GeF, only 
slightly above the D D* threshold. The experimental signatures for various 
ccud masses are also discussed. 
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I. INTRODUCTION 



Multiquark exotic hadrons different from tlie ordinary mesons or baryons liave been 
discussed and searclied for many years . Motivated by tlie possible discovery 

of doubly charmed ccu/ccd baryons, 0, we reconsider the cc m d tetraquark — the doubly 
charmed exotic state. We find that such a = 1^ state is likely to exist below or near the 
D* D thresholds, and may be the first stable/narrow exotic state to be discovered.0 

II. QQ'qq' STATES IN THE HEAVY QUARK LIMIT 

Let T{QQ' qq') denote a putative tetraquark state (it is denoted as Xc in [Q) consisting 
of two heavy quarks — Q, Q' = cot b — and two light quarks — q, q' = u, d. We are interested 
in genuine four quark "one bag states" with a "connected color network" . Figure |l] indicates 
one such state where Qa Q'p are combined via eap-y Qa Q'js into (3)^ color diquark, and q^ q'^ 
make an e^^^ q^ = (3)^ — an anti-diquark. The 3 and 3 combine to give an overall color 



singlet state. Note that according to the well-known arguments, |13;n|, in the large 
limit such a four-quark state is unstable against a decay into the two mesons, D and D* 
for instance. Alternatively, however, one can assume that the state T{ccud) in the large 
Nc limit corresponds to the state containing Nc — 1 heavy quarks, which combine into Nc 
representation of SU (Nc) color, and A^^^ — 1 light quarks combined into Nc. According to 
the standard large Nc counting rules the binding energy of such a {2Nc — 2)-exotic state, 
15| , [T6|1 , is of order 0{N^), i.e. of the same order as the meson masses. 



The identity e^/j-y e^r-y = ^aa ^pr — ^ar 5 (3a- allows one to express the above tetraquark 
state as a superposition of two meson states, {Q q) etc., which are separately color singlets. 

T{Q g' ^ 3, g r - 3} = I (Q q\ (Q' g')i ) - I (Q q')i {Q' q)i ) • (1) 



^Much effort was devoted in the past to search for hexaquark and pentaquark states [|,1C,11,12|. 
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FIG. 1. The color coupling pattern adopted here. Two charm quarks ci , C2 couple antisym- 
metricly to an intermediate 3 and u , d couple to 3. 



If the Q Q' (and q q') colors are coupled symmetrically to 6 (and 6) the analog of the Eq. (|I]) 
will have a plus relative sign. 

In a "string picture" the chromoelectric fluxes are squeezed into thin "vortices" connect- 
ing the various (anti)quarks and/or junction points. In this case the various lines in Fig. |l| 
describe not only the color coupling but also the actual layout of the strings. The tran- 
sition between the tetraquark (color connected — "one-bag" state) to the two- meson state 
can be pictorially described by shrinking the string bit connecting the two junctions, and 
then annihilating them via the above e ■ e contraction (Fig. |^). This naively would suggest 
the two-meson state, in which the above string bit has been eliminated, is lighter than the 
tetraquark state rendering the latter unstable since T{Q Q' qq') — > Q q + Q' q' would be 
kinematically allowed. 

However, the naive string picture may not apply to the ground state hadron considered 
here. Indeed we can directly show that T{Q Q' qq') is stable in the heavy quark limit 
{mq^Qi — s> oo). The Q and Q' would then bind into a 3 via the perturbative one gluon 
exchange. The essentially coulombic interaction yields a binding energy 0{a1)mQl2 (for 
iTiQ = mqr). Once ttlq is sufficiently large this binding exceeds hadronic energy scales and 
possible bindings in the heavy-light Q q mesonic systems. This then ensures the stability 
against decay into two such mesons. 

Unfortunately for this mechanism to generate stable tetraquarks rriQ (much) larger than 
(mc) nib is required. Detailed, and to some extent model dependent, considerations are thus 
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FIG. 2. The color flux vortices and their evolution as the tetraquark state separates into two 
D + D* states. 

required to motivate a. ccud state near or below D* D threshold. An alternative approach, 
which can directly address the stability of the physical 1"*" {ccud) state against a decay 
into D D* or D* D*, starts with these charmed mesons and attempts to form bound states 
via the potential due to light mesons, particularly one-pion exchange. Following Tornqvist, 
ir?! we will call such "deutron like" or "molecular" bound states of two mesons, "deusons" . 
The ranges and strengths of such potentials are independent of mq since D D* vr or D* D* ir 
couplings, for instance, depend only on the light quark degrees of freedom. Thus, if the two 
heavy mesons are attracted by these potentials, binding is again guaranteed in the heavy 
quark limit as the kinetic energy and centrifugal barriers for £ ^ waves vanish like I/uiq. 



III. DO D* D BOUND DEUSONS EXIST? 

Previous calculations utilizing OPEP (one-pion exchange potential) disagree on this is- 
sue. Thus Tornqvist [p!^ , |l^ , W\ finds B*B, D* D and even K* K bound states whereas 



Manohar and Wise |^ find only B* B bound states. The difficulty stems from the fact that 
much of the binding is due to the short range part of the potential.i The derivative pion 
coupling generates in the tensor part of the OPEP a e"^*"/ term which is singular at short 



^The central part of the OPEP in the D^*^ Z?*^*-* channel has the form: 

(ei • g) (e2 • q){ji ■ T2) 
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distance. Cutting off the OPEP at some distance ro is therefore required. The fact that 
Manohar and Wise (conservatively (?)) chose Tq ~ l/2mjr and Tornqvist (boldly (?)) takes 
ro ~ l/4m7r is the likely reason for their disagreement on bound D* D deusons. 

We do not believe that this issue can be settled. Yet the fact that the binding of this 
system by OPEP alone is not guaranteed may make the problem even more interesting. To 
decide the issue of a physical bound state with ccud flavor we may need the genuine four 
quark — "one bag" component of the state — where nontrivial aspects of QCD are operative. 
Indeed, once the distance between the D and D* mesons defined, say, by the distance between 
the respective charmed quarks, is smaller than the size of D or D* we can no longer treat the 
system as two separate hadrons exchanging light mesons. Rather, we need to revert to the 
"one bag" tetraquark description. The various quark- (anti) quark interactions in this state 
may then supply just the extra attractive interaction at (relatively) short distances needed 
in order to bind the system. 

The same physical c cm J state would then be a T(ccud) at short distances, r < ro(~ 
0.7 fm) and a. D* D deuson for r > Tq. Neither deusons nor tetraquarks are true ground 
states. The basic variational principle implies that mixing will lower the true ground state 
energy below the lowest energies found in each sector separately. 

IV. THE ccud TETRAQUARK 

Lacking a consistent first principle computational framework we appeal to the vast ex- 
isting lore and hterature. Thus to approach the problem of "color connected single bag" 
states one utilizes: 



with = — {ruD* — mo)'^ in D* — D and /x^ = in D* — D*. Averaging over the polarizations 
we get an expression like q^/iq^ + jJ^) ~ 1 — jic^ + yu^). The second piece becomes negligible as 
/X — >^ (except at g = 0) and the first piece contributes in configuration space a (5-function which 
clearly cannot be utilized in a reliable way for binding an extended object of interest. 
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i) "constituent" massive quarks ^ ^ 350 MeV, rric ~ 1.6 GeV ; 

ii) appropriate q/q~q long range interactions, or alternatively, an overall bag which confines 
the quarks into a single state; 

iii) the chromomagnetic hyperfine pairwise interactions: 

H.F. = n,, ^ — "|%(0)|2"(a. ■ (A. ■ A,) , (2) 

where cjj/Aj are the spin/color matrices of the (anti) quarks, and "|\&ij(0)p" is the relative 
wave function at zero separation for a qiqj meson and more generally the probability of 
overlap of the two (anti) quarks considered. 

Rather than attempting an ab initio calculation we adopt a more phenomenological 
approach. It utilizes known masses and insight from successful past calculations instead 
based on (i)-(iii) above in order to extrapolate to the ccud mass. 

We focus on the ccud I = state (rather then ccuu (ccdd) 1=1 states) since in 
both the tetraquark and deuson approach it is more strongly bound. The H.F. hyperfine 
interaction, Eq. (|^), strongly favors ud which is a 3 of color. The resulting anti-symmetry 
in color and separately in spin does then require fiavor antisymmetry, i.e. isosinglet u d state, 
in order to maintain overall Fermi statistics. Note that the ci C2 quarks, which are in 3 of 
color, must then form a spin triplet. This suggests an overall s-wave l"*" singlet ccud state 
whose lightest putative decay channel is indeed DD*. Also in the deuson approach the 
attractive OPEP is three times as strong in the / = than in the 1 = 1 channel. 

We thus consider the following putative double difference relation, 

("^(ccnJ)l+ - "^(cc«)l/2+) - ("^Ac - moo) =0. (3) 

It assumes that the extra energy required for replacing a u (or u) quark hj aud (ud) diquark 
in the presence of a c quark or (00)3 diquark is the same. It is clearly true in the heavy quark 
limit. Indeed in this limit the compact tightly bound Q Q pair is like a heavier antiquark 
flavor. Equation (^ would then be analogous to the relation inspired by the heavy quark 
symmetry: 
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(mA^ - tub) - {iriA, - ttld) = , (4) 

which holds very well. 

To further motivate Eq. we note that quark masses and almost all pairwise interac- 
tions cancel out via the double difference. The interactions between the Ci and C2 quarks 
in the cc = 3 diquark in ccud and in ccm match as well as the interactions between the 
ud = 3 (Is of spin) in ccud and the similar ud pair in cud. 

The H.F. interaction between the c quark and the u , d quarks in Ac = (cud) cancel out 
since (Tu = — (^d implies 

■ ^c) [K ■ Ac) + ■ (?c) (Ad ■ Ac) = ■ (Tc) ((A„ - Xdj ■ Ac) 
= -(a.-ac)(A^-A^) =0, (5) 

where we utilized the color neutrality condition, A„ + A^ + Ac = 0, for Ac. 

Similarly the H.F. interactions between the charmed quarks and the u , d quarks in 
T{ccud) cancel. Apart from an overall common factor \'$uc{0)\'^ /^ufnc we have, using 
again = - a^, 

i^ci ■ (?«) (Aci ■ Xu) + ((Tci ■ (Tj) (Aci • Aj) 

+ ('?c2 ■ <?«) (Ac2 ■ As) + {(fc2 ■ (?j) (Ac2 ■ Aj) 

= ((Tci ■ o^«) {Xu - Aj) ■ Aci + (o^c2 ■ (?«) (Xu - Aj) ■ Ac2 = Oa ■ (6) 

The last expression is an operator Oa anti-symmetric under the exchange of the color degrees 
of freedom of u and d. In the tetraquark state the u and d colors are also coupled anti- 
symmetrically in forming the 3 (anti) diquark. As a result the total H.F. interaction energy 
between the light and the heavy quarks in the tetraquark has the form (T4 \0a \ Ta ) which 
vanishes since upon an exchange of color indices of u and d this matrix element changes 
sign. 

To complete motivating Eq. (H) we still need to show that the cu (and cu) H.F. inter- 
actions in Sccu(l/2'^) and in D'^{0^) match. The latter is: 
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TUu rric 4 m„ 

where we used o^c + c?„ = and Ac + A^ = 0. 

The H.F. interactions in the doubly charmed baryon are however: 

i*«c(o)r ( 0"« ■ 0"cl) ( Au ■ Acl ) + (o"u ■ Crc2) ( Au ■ Ac2 ) ( 

1 l^uMl' . (^^^ ^ • (Aei + A,2)) 



2 m„ rric 

+ ■ (O^cl - 0'c2)) (A„ ■ (Xcl - Ac2)) } . (8) 

The second term above vanishes since, again, 

Au ■ (Aci — Ac2) ~ (^Aci + Ac2) ■ (Aci — Ac2) = — A^2 — 
so that we obtain: 

/ murric 2 murric 

where in the last step we used + 5^ + (?c2)^ = 3/4 (since the lightest charmed baryon 
has spin 1/2) and [ad + cTc2)^ = 2 (since the two charmed quarks are in state). 

We find that the hyperfine (attractive) interaction in Sccu(l/2^) and in D°(0~) do not 
exactly match but are slightly stronger in by 

^ ''H.F." = J {mD" - niD) ■ 

Thus by subtracting the physical lighter D mass we are actually causing an imbalance in 
Eq. (^ and an overestimate of rrirpf^^^^^y A corrected equation should therefore read: 

rriTiccud) = ms^cu + - - ^ {rriD' - mn) ■ (10) 

While the masses of Ac, and D^* are well known, this is definitely not the case for the 
doubly charmed baryons. The lowest SELEX peak appears in ccu*-"*"^-' at 3460 MeK 0. 
Using this value in Eq. (|iy) yields 

^T(cc^J)~ 3.845 Ge\/, (11) 



which is about 25 MeV below the D* D threshold. 

Unfortunately the lowest peak in ccS^^ is at 3.52 GeV and not degenerate with the ecu 
peak as it should be by isospin invariance. We therefore choose the ccS~^^ peak (which 
indeed coincides with a (relatively small) enhancement in the ccrf*-^-* SELEX data) as rep- 
resenting the true value of the lowest charmed baryon. This is clearly a more conservative 
choice as both SELEX peaks are lower than most previous theoretical predictions. Using 
this we obtain, 

^T(cc^<J) - 3.905^61^ (12) 

some 35 MeV above the D* D threshold. 

Various explicit and implicit assumptions were made in order to obtain Eq. ( p!OD mak- 
ing for a theoretical uncertainty in the predicted value of my(cctij) beyond the 60 MeV 
experimental uncertainty discussed above. 

a) We used a common overlap probability "|\E'jj(0)p" for all eq or eq pairs which cor- 
responds to a "universal bag radius" . In reality the latter could change as we go from two 
to three to four quark systems. The successful phenomenology of baryon/meson hyperfine 
splitting suggests that this may be a weak effect. Furthermore a systematic change will 
largely cancel in the double difference relation in Eq. (|^). 

b) We have restricted our discussion to diquarks coupled to 3 (or 3) of color only. However 
in the tetraquark system we encounter (for the first time) the possibility of coupling Ci C2 
(and ud) to 6 (6) of color — a coupling, which is clearly not allowed in a baryon. The 
restriction to the 3 — 3 pattern may be well justified by the fact that it has a lower energy 
than the 6 — 6. Nonetheless these channels can in principle mix, and allowing for such 
admixture to optimize the binding will lower the mass of the physical state. 

c) The admixture of the "one bag" tetraquark state above and the "two bag" deuson 
state (which again occurs first in four-quark exotic states) will also, by the above mentioned 
variational argument, tend to lower the energy. 

All the above suggests that if the SELEX peak at 3460 MeV is indeed the lightest 
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double charmed baryon, we have a T(c cud) state shghtly below or slightly above threshold. 
Specifically, if 

(^ = mT(ccud)-mD* -tud, (13) 

we expect 

|e| < 30 ^60 Mel/. 

V. PRODUCTION AND DECAYS OF T{ccud) 

We would next like to argue that if T{ccud) has the above mass it may well be the first 
narrow exotic hadron to be discovered. The potential discovery depends jointly on 

(i) the rate of c c -u d production [|22|, and 

(ii) the existence of decay modes which can provide a unique signature. 

The production rates (at hadronic or colliders) of the state ccud of interest are 

very small as all the following conditions should be met: 

(a) Two pairs of charmed quarks (ci ci) and (c2 C2) need to be produced. 

(b) These pairs should be close spatially. Also Ci from the first pair, and C2 from the second 
should have small relative momenta in order for a a Ci C2 diquark to form. 

(d) Finally the Ci C2 diquark should pick np a ud (anti)diquark to form ccud. 

The first two factors (a) and (b) also suppress the production rate of doubly charmed 
baryons ccu/ccd. The only further suppression of the ccud production rate is due to the 
need to pick up a diquark instead of merely just one u (or d) in the case of ccu/ccd. 
This suggests the following "double ratio" relation: 

{R{ccud)/ R{ccu)) ■ [R{cud)/R{cu)) = 1, (14) 

which is analogous to the double difference relation in Eq. (|^).i 
^This analogy is very natural in a "statistical model" , where the production rate of any hadronic 
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The ratio of the charmed baryon, production and that of D's is roughly the same as 
in the case of strange baryons/mesons: 

R{cud)/R{cu) ~ R{sud)/R{su) = R{K)/R(K) ~ Yq • 



Hence we expect from Eq. (0) 



^{ccud) ~ Y^Riccu) , (16) 



namely that the T{ccud) production rate is about 1/10 that of charmed baryons. 

If T{ccud) is to be discovered this tiny production rate needs to be compensated by 
striking decay signatures. The decay modes critically depend on rrirpf^^^^^-^ or e, the separation 



between rrirp^^^^j^^^ and the rriD + rrij:,* threshold, defined in Eq. (|1^). 

We will next discuss the different mrpi^^^^g^ ranges starting with the most strongly bound 
case. 

(a) J) < 2 m^) oi e < m£, — mjj* ^ — 140 MeV. 

In this case — which the above discussion suggests to be unlikely — we can only have two 
consecutive weak decays. Strictly speaking only the second vertex involves an on-shell 
reconstructible or D^. However, up to small binding effects, 2 itid — rrirp^^^^^g^, the tracks 
emerging from the first decay may well correspond to another D. 

There is some probability that the charm quark surviving after the first decay may be 
inside a D* so that an extra, slow pion emitted from the first decay vertex can combine with 
the 4-momenta of the second vertex particles to form a D* . Finally all particles from the 
two decay vertices should reconstruct a narrow T{ccud) state below 2mD. 

(b) rriD* + rriD > mT(^ccud) > 2 ttid or > e > — 140 MeV 
In this case we will have an electromagnetic decay: 

T(ccMrf)+ ^D+ + D° + 7, (17) 



state X is suppressed by a Boltzmann factor proportional to exp {—mx/T). Equation (p!^) results 
then as the exponential of Eq. (P). 
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with two weak decay vertices due to D'^ and D° decays. While the state will be narrow, 
reconstructing the invariant mass, 171^+ jjOy, and looking for a sharp peak requires identifying 
a relatively soft photon [E^ < 140 MeV in the T[ccud) rest frame) emerging from the 
primary decay vertex. 

(c) rriTf^^^r^S) > TTT'D* + or e > 
Here we clearly have T{ccud) D*^ + D+ (or D*+ + D^). The D* will decay into D + tt 
at the primary vertex and the two D's will next decay weakly at separate vertices. Thus 
we attempt to reconstruct the two D^s and the D* and finally look for an overall peak in 
the D* D invariant mass distribution. This last constraint may not be very helpful (and 
the T{ccud) can be altogether missed) if its width r{T{ccud)) (substantially) exceeds the 
experimental resolution which we optimistically take to be O{10 MeV). 

Since T{ccud, 1"*") D* D is an ^-wave decay and no new quark pairs need to be created 
(as in K* — > Kn, for instance) one might expect a large decay width, F 300 MeV, as 
is the case with qj qi exotics made of light quarks. Two factors may, however, reduce 
r{T{ccud)). First we have a two-body decay phase space which is proportional to /3 or 
the velocity of D or D* in the T{ccudys rest frame, which unhke in the decays of light 
exotics may be significantly less than one: 

Second, it may well be that the physical l"*" hadron of interest has a relatively small 
deuson component, \a\'^: 

\^ ^^ji^ (physical) ) = \Jl — jap \T{ccud)) + a \ {D* D) deuson) , (19) 

with T{ccud) being the "one bag" genuine four quark state. Since only the deuson com- 
ponent readily falls apart into D* + D we may then have a further suppression by a factor 
of \a\^. If \a\^ < I — \ the joint \a\^ (3 effect reduces the decay rate from F fa 300 MeV to 
Ffa 13 -19 MeV. 

Note, since the state of interest is an S'-wave of the D + D* we may lack any repulsive 
interaction — akin to Coulomb repulsion in fission or a decays — ^to generate a resonance in the 
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first place [^. Rather we may have, as in the 1 = 1 S-wave nucleon-nucleon scattering, a 
"virtual bound state" . The experimental manifestation of the latter — a strong enhancement 
at the D D* threshold, may be still sufficiently striking. The above lifetime estimates would 
then still pertain to the width of this enhancement. 

In principle we can also have cued tetraquark (and/or meson-antimeson, D* D, D D* 
deusons) with "hidden" charm. Such combinations will be much more easier to produce 
as only one cc pair needs to be created. Naively one would expect such states to decay 
very quickly into J/\E' + vr^. Tornqvist, whose main concern was actually meson-antimeson 
binding did, however, suggest that a primary extended deuson D* D state would have little 
overlap with the J/\l/ (compact cc) state. This would then make for a relatively narrow 
D* D state which still could be nicely identified via its unique {J + vr+) decay mode at a 
specific J/\l' + 7r+ invariant mass. 

The ccu/ccd double charmed baryons have presumably been seen at FNAL 0. The 
T{ccud) is most likely produced there with high lab momenta. The proximity of the 
T{ccud) mass to the D* D threshold cause the two subsequent weak decay vertices to be 
very precisely aligned with the initial interaction vertex — a feature which will be most helpful 
in a tetraquark hunt. 

The much lower energy e~ colliders can also serve as promising search grounds with a 
much cleaner environment. Thus at ^/s ^ 10.7 GeV (~ wix(4s)) the primary virtual photon 
interaction yields in about 25% of all cases a cc pair. This then will vastly increase the 
number of events in which we have two charmed quark pairs to be 0{1 — 0.25), i.e. O{10^) 



events [^. The events in which (J/\l/) will not form will have typical 



final states (D^*^ indicates D or D*) with a few extra particles due to the limited phase 
space (4m£)* ~ 8.4 Gel^!). The systematic search of peaks in D* D^*^ or D^*^ D invariant 
mass distribution could thus be feasible despite the large combinatorial background. 

Peaks in D* D^*^ mass distribution could also occur. These should manifest in the much 
cleaner (J/\E' + vr) channel. Hopefully by D* (or D**) cascade decays and by K^, K~ 
separation one will be able to distinguish D^*^ D^*^ from D^*^ D^*^ pairs. 
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Finally we note that the Brown- Hanbury-Twiss effect favoring D° (and D'^) pairs 
with small relative momenta should not be operative here as the two pseudoscalar D's do 
not emerge from the same vertex (even the D's from the D* ^ Dtt decay would emerge 
several hundred Fermies away from the primary vertex because of the narrow D*'s width). 

VI. SUMMARY AND CONCLUSIONS 

We have presented various estimates pertaining to a possible tetraquark T{ccud){l~^) 
state. If the SELEX second peak corresponds to the lightest doubly charmed baryon, then 
our estimates of the mass of T{ccud){l~^) are close to the D D* threshold. We have also 
discussed signatures for various possible masses. 
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